home *** CD-ROM | disk | FTP | other *** search
/ Language/OS - Multiplatform Resource Library / LANGUAGE OS.iso / cpp_libs / rwvector.lha / RWVector2.1 / src / mathpack / cffti.f < prev    next >
Encoding:
Text File  |  1989-08-14  |  1.4 KB  |  46 lines
  1. *deck cffti
  2.       subroutine cffti (n,wsave)
  3. C***BEGIN PROLOGUE  CFFTI
  4. C***REVISION DATE  811015   (YYMMDD)
  5. C***CATEGORY NO.  D6
  6. C***KEYWORDS FFT,FAST FOURIER TRANSFORM,FOURIER TRANSFORM,COMPLEX
  7. C***DATE WRITTEN  FEBRUARY 1978
  8. C***AUTHOR  SWARZTRAUBER P.N. (NCAR)
  9. C***PURPOSE
  10. C   INITIALIZE FOR CFFTF AND CFFTB
  11. C***DESCRIPTION
  12. C     *****************************************************************
  13. C
  14. C     SUBROUTINE CFFTI(N,WSAVE)
  15. C
  16. C     *****************************************************************
  17. C
  18. C     SUBROUTINE CFFTI INITIALIZES THE ARRAY WSAVE WHICH IS USED IN
  19. C     BOTH CFFTF AND CFFTB. THE PRIME FACTORIZATION OF N TOGETHER WITH
  20. C     A TABULATION OF THE TRIGONOMETRIC FUNCTIONS ARE COMPUTED AND
  21. C     STORED IN WSAVE.
  22. C
  23. C     INPUT PARAMETER
  24. C
  25. C     N       THE LENGTH OF THE SEQUENCE TO BE TRANSFORMED
  26. C
  27. C     OUTPUT PARAMETER
  28. C
  29. C     WSAVE   A WORK ARRAY WHICH MUST BE DIMENSIONED AT LEAST 4*N+15
  30. C             THE SAME WORK ARRAY CAN BE USED FOR BOTH CFFTF AND CFFTB
  31. C             AS LONG AS N REMAINS UNCHANGED. DIFFERENT WSAVE ARRAYS
  32. C             ARE REQUIRED FOR DIFFERENT VALUES OF N. THE CONTENTS OF
  33. C             WSAVE MUST NOT BE CHANGED BETWEEN CALLS OF CFFTF OR CFFTB
  34. C
  35. C***REFERENCES
  36. C***ROUTINES CALLED  CFFTI1
  37. C***END PROLOGUE  CFFTI
  38.       dimension       wsave(1)
  39. C***FIRST EXECUTABLE STATEMENT  CFFTI
  40.       if (n .eq. 1) return
  41.       iw1 = n+n+1
  42.       iw2 = iw1+n+n
  43.       call cffti1 (n,wsave(iw1),wsave(iw2))
  44.       return
  45.       end
  46.